Giant Surface Plasmon Induced Drag Effect (SPIDEr) in Metal Nanowires 
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Here, for the first time we predict a giant surface plasmon-induced drag effect (SPIDEr), which 
exists under conditions of the extreme nanoplasmonic confinement. Under realistic conditions, in 
nanowires, this giant SPIDEr generates rectified THz potential differences up to 10 V and extremely 
strong electric fields up to ~ 1 5 — 1 6 V/cm. The SPIDEr is an ultrafast effect whose bandwidth for 
nanometric wires is ~ 20 THz. The giant SPIDEr opens up a new field of ultraintense THz nanoop- 
tics with wide potential applications in nanotechnology and nanoscience, including microelectronics, 
nanoplasmonics, and biomedicine. 



i 

C 

o 
o 



(N 

o 

OS 

o 



X 



The fact that electromagnetic fields exert mechanical 
forces on matter is well known and has found applications 
in atomic physics-^ and bio- and nanotechnology; 3 -^ 
and picosecond photodetectors based on the photon drag 
effect (PDE)i£ The semiconductor PDE detectors have 
proved to be very practical for relatively fast detection 
of picosecond pulses in a wide frequency range spanning 
from THz to infrared. It has beeen proposed to use en- 
hanced fields in phonon-polariton silicon carbide struc- 
tures for laser particle accelerators^ There have been ex- 
perimental investigations of the PDE in metal nanofilms 
that are thicker than the skin depths Metals support 
surface plasmon polaritons (SPPs)i2. that exert forces on 
electrons causing an SPP-enhanced PDE.» However, it 
is typically rather small, with the induced potential dif- 
ferences in the mV range. This modest enhancement 
of the PDE in these experiments was due to relatively 
slow variation of SPP fields in space, and the momentum 
transferred to the electrons was correspondingly small. 
Because of the long relaxation times of the SPPs at plane 
metal surfaces, this PDE is also expected to be rela- 
tively slow. A drag effect under the conditions of strong 
nanoplasmonic confinement, when the SPP localization 
radius is less than the skin depth (~ 25 nm), has not been 
studied or exploited theoretically or experimentally. 

In this Letter, we predict a giant surface plasmon- 
induced drag effect (SPIDEr) in metal nanowires, which 
is fast, with response on the femtosecond time scale. We 
show that the ultrashort, nanolocalizcd SPP pulses exert 
forces on electrons in the nanowires, inducing giant THz 
electromotive force (emf) along the SPP propagation di- 
rection. We have found that in thin (~ 5 nm radius) 
wires this emf can reach ~ 10 V, with nanolocalized THz 
fields as high as ~ 1 MV/cm. Such THz field have previ- 
ously been generated in the far zone 1 !., where they pro- 
duce non-perturbative effects^ but not on the nanoscale. 
In contrast, the plasmonic metal nanowires can serve as 
nanolocalized sources of high THz fields. We also study 
the dynamics of the SPIDEr for short SPP pulses and 
suggest that adiabatically tapered nanowires^ can be 
used as broadband nanoscopic photodetectors with ex- 
tremely fast response due to the femtosecond momentum 
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FIG. 1: Schematic of SPIDEr in metal nanowire. Propagating 
SPPs create forces acting on carriers in the nanowire, which 
leads to THz-band voltage (emf) between the ends of the wire. 
Picosecond or femtosecond pulses can be used to manipulate 
the time dependence of the created emf. 



relaxation times in metals^ The nature of the giant en- 
hancement of the SPIDEr is novel in nanoplasmonics: it 
is not the enhancement of the optical fields per se (the 
maximum magnitude of the local fields is limited by the 
breakdown at the metal surface that occurs for fields ~ 1 
V/A) but the very high gradients of these fields. The 
SPIDEr is ultrafast because it is a non-resonant effect 
whose bandwidth is comparable to that of the entire op- 
tical spectrum. 

Consider a metal nanostructure propagating an SPP 
pulse. The SPP field produces polarization P, charges 
with macroscopic density p = — divP, current density 
j = dP/dt, and the surface charge density a — (P • n) 
at the surface of the metal, where n is the normal to the 
surface pointing outward. We do not consider systems 
with optical magnetism, whose introduction at optical 
frequency is problemati o 15 i 16 Therefore we set B = H, 
which precludes the existence of surface currents. Us- 
ing Eqs. ©-© of the Methods Section, we obtain the 
following general expression for the total force 
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where superscript "c" implies that the differentiation 
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does not apply to the labeled vector. This result is of fun- 
damental importance for processes involving interaction 
of nanoplasmonic fields with metal electrons.— Equation 
(fTj) is valid for a wide range of problems with a general 
material equation P = %E, including those where opera- 
tor x describes anisotropic or non-local media. The first 
term in Eq. ([1]) is similar to the force acting on a point 
dipole moment.-- The second term in Eq. (JT|) is the Abra- 
ham force. In a monochromatic field, this force averaged 
over the period of oscillations is zero, but in the field of 
a pulsed excitation it has a finite magnitude. 

We apply this fundamental Eq. |T]) to describe the SPI- 
DEr. For certainty, consider a metal nanowire with ra- 
dius R and dielectric susceptibility x( w )j which is ori- 
ented along the z-axis and embedded into a dielectric 
with a dielectric permittivity of Ed- This wire propagates 
an SPP pulse, which can be excited by external sources 
using, e.g., the effect of adiabatic compression . 13 ! 18 

In the case of extreme nanoplasmonic confinement 
(R <C i s , where l s is the skin depth), R becomes the 
only relevant quantity of the dimensionality of length. 13 
Therefore there is scaling of all magnitudes in R. The 
SPP wave power £P scales as ~ E 2 R 2 v g , where v g ~ 
loR is the SPP group velocity. The SPIDEr-induced po- 
tential difference [electromotive force (emf )] S is propor- 
tional to the pressure produced by force fl}, S ~ F/R 2 . 
The propagation length of the SPP l p ~ RQ, where Q 
is the SPP figure of merit, independent of R. These ar- 
guments allow us to predict scaling of the SPIDEr force 
F, emf S ', and the electric field due to SPIDEr rectifi- 
cation Eji (which for femtosecond SPP pulses possesses 
THz frequencies): F oc ^R' 1 , S oc ^R~ 3 , E R oc 
£PR~ A , E m R oc i?" 1 .We have also indicated the scal- 
ing of the maximum rectified field E m R (at the maxi- 
mum tolerable power &m). The scaling implies that all 
the effects caused by the SPIDEr increase with decreas- 
ing the wire radius as its powers. This enhancement is 
not resonant and therefore has bandwidth comparable to 
that of the entire optical spectrum. The scaling describes 
only the dependence on R. There are also prefactors de- 
scribing the dependence on dielectric permittivities, fre- 
quency, etc. They take into account an additional en- 
hancement close to the SP resonant frequency, which is 
multiplicative in magnitude. Below in this paper we show 
that this scaling is reproduced by the theory results. 

The SPPs are transverse magnetic (TM) modes, and 
their complex fields have the form 

E = A(i!)(E z i+E p py {kz -" t] , H = Ait^H^e 1 ^-^, 

(2) 

where t = t — z/v g , and v g is the SPP group veloc- 
ity at the pulse carrier frequency u>, k is the SPP wave 
numberJ^ The total power flowing through the plane 
z = at the moment t is &{t) — 2ir J °° pdpS z (r, t)\ z =o, 
where S(r,t) = (c/87r)Rc[E x H*] is the Pointing vector 
averaged over SPP period. Considering the azimuthally- 
symmetric (TM ) modes, functions E z , E p and H v de- 
pend only on radius p. We normalize them for real ampli- 
tude A(t) to satisfy a relation A 2 (t) — 3P(t). Disregard- 



ing the group velocity dispersion in Eq. ([2]) is valid for 
pulses with duration of tens of femtoseconds and greater, 
and frequencies not too close to the SP resonance. 

The momentum transferred from the electromagnetic 
field to the electronic system implies forces exerted on 
the electrons. The density of these forces are given by 
Eqs. and {5| of the Methods Section. These forces 
lead to an emf, which corresponds to optical rectifica- 
tion in the system, which for femtosecond SPP pulses 
results in the emf in the THz frequency range. Since 
the electron momentum-relaxation time is on the scale 
of femtoseconds, electrons come to a local equilibrium in 
the process of this rectification. Therefore, we will de- 
scribe it in the hydrodynamic approximation, for which 
the pressure p and electrostatic potential <j> satisfy an 
equation p + ne<j) — const, where e is electron charge, 
and n is electron density. From this equation, we can 
find the emf S — A(f>, which is the total change of poten- 
tial in the direction of SPP propagation (the z direction), 
<8 = —Ap/(ne), where Ap = F z /(irR 2 ) is the full change 
of the pressure. Here, F z is z-component of the force ([T]) 
averaged over the period of SPP oscillations. 

The total force F z is composed of three forces (see Eq. 
(|9]) in Methods Section): the SPP pressure, striction, and 
Abraham force F z = ff r + jf + f A . These three forces 
result in three terms of the SPIDEr emf 

S R n I m I *'®) , (3) 

V A pr A st cL A J 

where Rh = — l/(ecn) is the Hall constant, and coeffi- 
cients Ap r , A s t, and La and &(t) arc defined by Eqs. 
pQ ]) -(fT3 ]) of the Methods Section. The Abraham force 
contribution [the third term in Eq. (|3|)] is small under 
a condition kolp/{kcr) <C 1, where fco = w /c, and r 
is the pulse duration. This condition is well satisfied 
for the parameters used in this paper. We will con- 
sider two limiting cases pertaining to Eq. ©: a quasi- 
monochromatic regime of long pulses (r ^> t p ) and a 
regime of short pulses (r <C t p ), where t p is the SPP dis- 
sipation time. Consider first the quasi-monochromatic 
regime where &(t) « 3 s {t) (see Eq. (TTU]) of the Meth- 
ods Section). In such a case, the total emf $(t) follows 
the pulse-envelope time dependence 3 s (t). The compu- 
tations will be made for a silver— nanowire embedded in 
vacuum (e^ = 1). 

For the quasi-monochromatic case, the total emf § as a 
function of frequency to and wire radius R is displayed in 
Fig. [5] (a) . In contrast to the case of dielectric media, in 
the present plasmonic case the SPP pressure and stric- 
tion contributions to the emf have the same sign since 
x' < 0. These two contributions are equal if a condition 
x"Q = —x' is satisfied. The black solid line represents 
this condition; to the left of this line the pressure domi- 
nates, and to the right the striction force gives the major 
contribution to the SPIDEr. This is understandable be- 
cause close to the SP resonance of the wire (at w 3.7 
eV), the gradient of the SPP intensity increases due to 
the high loss: the striction force is of a gradient nature, 
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FIG. 2: SPIDEr for quasi-monochromatic SPP pulses: emf 
and rectified field dependence on the frequency hut and wire 
radius R. Note the logarithmic scale for the magtnitude of 
the effect, (a) Dependence of the SPIDEr emf per unit SPP 
power £ I on wire radius and frequency. The black broken 
curve indicates the parameters at which SPP pressure is equal 
to striction. The magnitude of the effect is denoted by the 
color-coding bar. (b) Dependence of SPIDEr magnitude S 
on wire radius R per unit power of the SPP wave (solid red 
curve). The contributions of the pressure and striction to the 
total magnitude of SPIDEr are shown by the dashed curves, 
(c) Dependence of SPIDEr magnitude per unit power $ / '8? 
on frequency w for R = 5 nm. (d) Maximum power that a 
wire can tolerate S?m as a function of wire radius R for differ- 
ent frequencies u>. (e) The maximum SPIDEr rectified field 
E m R (for the maximum tolerable power & m ) as a function 
of frequency for three wire radii R = 5, 30, and 100 nm. (f) 
The maximum SPIDEr rectified field E m R (for the maximum 
tolerable power as a function of the wire radius R for 

frequency hiv — 1.55 eV. 



dominating therefore. Similarly, with the decrease of R, 
the intensity gradient increases due to the increased con- 
finement, which also leads to the relative increase of the 
striction with respect to the pressure force, as we clearly 
can see from this and other panels of Fig. O General 
increase of the SPIDEr at the SP resonant frequency can 
be seen as a broad red peak. 

Magnitude of the SPIDEr emf relative to the SPP wave 
power, S '/ ' 3? ', is illustrated in Fig. [3] (b) as a function 



of the wire radius R for a frequency of fuo = 2.9 eV. 
The SPIDER effect is gigantically enhanced for strong 
nanoplasmonic confinement: by four orders of magnitude 
when R decreases from 100 to 5 nm at the same SPP 
power. There is a pronounced scaling <S j 3? oc R~ 3 at 
R< l s , in accord with the discussion following Eq. (JTJ) . 

The spectral dependence of the relative SPIDER emf, 
£ ' / '3 d ', for a wire of the smallest radius considered, R = 5 
nm, is depicted in Fig. [He). 27 Importantly, the magni- 
tude of the emf in this case is very large, from 0.01 to 10 
V per 1 W of the SPP pulse power, in the entire optical 
range, with a pronounced resonance at the SP frequency. 
This large magnitude shows that the SPIDEr effect can 
be used for the photodetection on the nanoscale, i.e., in 
the role that previously was deemed only available for 
semiconductors. In this sense, it belongs to the area of 
the active nanoplasmonics.- 20 

By classification of the nonlinear optics, the SPIDEr is 
an optical second-order nonlinear effect: the magnitude 
of the SPIDEr emf & is proportional to the power 3 d of 
the SPP field. Therefore, the maximum achievable mag- 
nitude of the emf is determined by the maximum 3 that 
the wire can tolerate. This we estimate setting the opti- 
cal field E at the surface of the wire equal to 1 V/Ai 21 i 22 
For fields significantly higher than this, there will be mas- 
sive ionization and damage of the metal surface. We plot 
in Fig. [2] (d) this maximum intensity as a function of the 
wire radius for three SPP frequencies. Note a very good 
scaling ^ m oc R 3 in the region R < l s of the strong 
nanoplasmonic confinement. These values of the 3 s m in 
comparison to the data of Fig. [5] (a)-(c) show that the 
gigantic values of of the SPIDEr emf S ~ 10 V are real- 
istically achievable, which are many orders of magnitude 
greater than observed previously in the metal filmsi^ 

One of the most important for applications properties 
of the giant SPIDEr is high local electric field En gener- 
ated due to the SPIDEr optical rectification in the vicin- 
ity of the nanowire. Such a field (averaged over the SPP 
decay length l p ) can be found as En = S jl v . We dis- 
play the maximum achievable rectified field E m n (at the 
propagating SPP power of 3 s m ) in Figs. [5] (e), (f). As we 
see from panel (e), the spectral dependence of the SPI- 
DEr rectified THz field is very similar for all wire sizes, 
but the magnitude of this field is much greater for the 
5 nm wire: E m R ~ 10 5 — 10 6 V/cm. The nanolocalized 
THz fields of such a magnitude will excite a wealth of 
nonlinear THz responses at the nanoscale. 

Now let us consider dynamics of the emf response to 
SPP pulses that differ in duration r with respect to the 
the SPP pulse dissipation time t p . The latter is displayed 
in Fig. [3] (a) as a function of the frequency u>. As we 
can see time t p is in the range from 10 to 150 fs. The 
temporal dependencies of the emf in comparison with the 
power 3 d of the SPP pulses for various pulse durations is 
illustrated in Figs. [3] (b)-(d). For a relatively long pulse 
(t = 1 ps > t p ) shown in Fig. [3] (b), the shape of the 
emf (o(t) repeats that of the power 3 s it). This relatively 
long, picosecond response, nevertheless, corresponds to a 
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FIG. 3: SPIDEr created by ultrashort SPP pulses: fast fem- 
tosecond emf response, (a) The dependence of the SPP life- 
time t p = lp/vg on the frequency hu for R — 5 nm and 
R — 30 nm. (b) The time dependence of the emf £{i) (green 
line, left scale) and input power ^(t) (red line, right scale). 
The pulse duration is r = 1 ps S> t p ~ 30 fs and the emf 
closely follows the SPP pulse dynamics, (c) The same for 
much shorter pulse with r = 40 fs. The pressure-induced 
emf leads to small broadening in emf dynamics (green line). 
The limiting exponential decay is outlined by the broken blue 
line, (d) Emf induced by the short pulse in nanowire with 
R = 30 nm with frequency huj = 1.2 eV. The emf response is 
broadened, since r = 80 fs, while t p ~ 150 fs. 



1 THz bandwidth for this nanowire used as a nanoscale 
photodetector. Note that the amplitude of the emf is 
very large, ~ 10 V. 

For a much shorter, r = 40 fs, SPP pulse and the same 
5 nm wire, as shown in Fig.[3](c), there is a small broaden- 
ing and delay of the voltaic response (emf) <§ (t) with re- 
spect to the excitation SPP pulse ^(t). This broadening 
is due to the pressure force that decays exponentially for 
long times, as Eq. (fT0|) of the Methods Section suggests, 
and the broken blue line in the figure indicates. However, 
under the conditions considered, this delay and broaden- 
ing are not large. The frequency-response bandwidth of 
this wire as an SPP photodetector on the nanoscale is 
very large, « 20 THz, which is characteristic of the ex- 
treme nanoplasmonic confinement. The amplitude of the 
SPIDEr emf is also very large, S ~ 10 V. 

For a much thicker nanowire of R = 30 nm (weak 
plasmonic confinement case) and fku = 1.2 eV, illus- 
trated in Fig. [3] (d), the SPP decay time becomes much 
longer (t p = 150 fs). This leads to a very significant de- 
lay and temporal broadening of the emf response with a 
pronounced exponential part due to the pressure forces 
shown by the broken blue line. These behavior is due to 
the much longer SPP lifetimes for the weak confinement 
where a significant fraction of the SPP energy propagates 
in the dielectric (vacuum). Nevertheless the emf response 



bandwidth is still very large, on the order of 5 THz, and 
its amplitude is also very large, S > 1 V. 

In conclusion, the ultrafast giant SPIDEr in metal 
nanowires excited by ultrashort SPP pulses is predicted 
in this Letter to generate a gigantic emf up to ~ 10 V 
for the SPP waves of realistic and tolerable amplitudes. 
The SPIDEr enhancement is mostly nonresonant, due 
to nanoplasmonic confinement of the SPP fields, which 
leads to higher gradients of the fields and their higher 
magnitudes. Because of its nonresonant nature, the SPI- 
DEr is an extremely fast effect: frequency bandwidth of 
the generated THz fields is realistically 5-20 THz. Due 
to high longitudinal localization of the SPP waves in the 
case of the strong nanoplasmonic confinement, the SPI- 
DEr generates very high local THz electric fields at the 
metal surface, Er ~ 10 5 — 10 6 V/cm. Such fields are ca- 
pable of inducing strongly nonlinear responses, including 
dissociation of molecules. Among possible applications 
of the giant SPIDEr are rectification and detection of 
the nanoscale femtosecond optical fields, coupling of the 
nanoplasmonic elements to semiconductor devices, non- 
linear THz spectroscopy on the nanoscale of chemical and 
biological nanoobjects for biomedicine, etc. 
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Methods 

Then a macroscopic (averaged) Lorentz force density 
inside the material volume is given by 

f = -E(V-P) + -JxB. (4) 

c at 

The surface polarization charges experience the Lorentz 
force with a surface density 

f 1 = crE . (5) 

The volume force density given by Eq. (U) can be repre- 
sented in the following way 

r = -div (P ® E) + grad(P c ■ E) + - ^ ' * B) , (6) 

c at 

where ® denotes the outer product of vectors. Thus, the 
total force acting on the polarization charges in the SPP 
field is 

F = / f" dV + <£ f s ds , (7) 
Jv Js 
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where the first integration runs over volume V of the 
structure, and the second is over its surface S. Noticing 
that the volume integral of the first term in Eq.® is 
equal to 



/ div (P ® E) dV = - * (P ■ n)E ds = - I f s 
Jv Js Js 



ds 



(8) 

we cancel terms in Eq. to simplify it and obtain the 
required Eq. (1) of the paper. 

From Eq. (1) of the paper, the total force acting on 
electrons can be expressed as a sum of the SPP pressure, 
striction, and Abraham force, whose ^-components are, 
correspondingly, 



nR 2 J?(t) 



cA 



fP = 



> -nR 2 &'{t) 



pr 



cA, 



c 2 L A 



(9) 



where 



1 

lp Jo 

e -t/t P 

trj 



^(t')e- z/l "dz = 
^{t')e t,,tp dt' . 



In Eq. |9]), coefficients A pr and A st have dimensionality 
of area, and La has dimensionality of length. They are 
defined as 



La= \%Iv Re {*E P H;)pd P 



(ii) 

(12) 
(13) 



where Q = Refc/Imfc is the SPP figure of merit, |E| 



\E Z \ 2 , l p = l/(2Imfc) is the SPP propagation 



(10) length, and t p — l p /v g is the SPP pulse lifetime. 
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